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ABSTRACT 



We have fitted ~ 200 RXTE and INTEGRAL spectra of the neutron star LMXB 
GX 9+9 from 2002-2007 with a model consisting of a disc blackbody and another 
blackbody representing the spreading layer (SL), i.e. an extended accretion zone on 
the NS surface as opposed to the more traditional disc-like boundary layer. Contrary 
to theory, the SL temperature was seen to increase towards low SL luminosities, while 
the approximate angular extent had a nearly linear luminosity dependency. Comp- 
tonization was not required to adequately fit these spectra. Together with the ~ 70° 
upper bound of inclination implied by the lack of eclipses, the best-fitting normaliza- 
tion of the accretion disc blackbody component implies a distance of ~ 10 kpc, instead 
of the usually quoted 5 kpc. 

Key words: stars: individual: GX 9+9 - X-rays: binaries - accretion, accretion discs. 



1 INTRODUCTION 

The Galactic X-ray source GX 9+9 (3A 1728-169, 4U 1728- 
16) was discovered by a sounding rocket fligh t perfo rmed 
on 1967 July 7, as reported by iBradt et all l|l968f ). An 
optical counterpart of magnitude V = 16.6 with a blue 
continuum and a larg e ultraviolet excess was identified by 
iDavidsen et al.l ljl976l ). 

GX 9+9 was found to have similar spectral and tim- 
ing properties to Sco X-l and other sources that came 
to be called Low-Mass X-ray Binaries ijMason et al.l Il976l : 
IParsignault fc Grindlavlll978l ). It was classified as an atoll- 
ty pe neutron star source in a pers istent lower banana state 
bv lHasinger fc van der Klid ljl989h . 

The discovery of an X-ray modulation period from 
observations performed by the High Energy Astronomy 
Observatory- 1 (HEAO -1) in 1977 September was announced 
llHertz fc Woodl Il986h and refined to 4.19+0.02 h by 
iHertz fc Wood! l|l988h . The p eriod wa s inter p reted as the 
binary period of the system. ISchaeferl ljl987l . Il990h found 
a similar period of 4.198+0.0094 h in the optical counter- 
part. This was attributed to the orbital modulation of the 
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X-ray emission reprocessed into visible light in the compan- 
ion star atmosphere and at th e hot spot wher e the a ccretion 
st ream hits the ac cretion disc. lHertz fc Woodl l|l988l ). as well 
as ISchaeferl fl990b deduced the companion to be an early M- 
class dwarf, with mas s estimates of 0.2- 0.45 and ~0.4 Mq, 
respectively. However, iReig et all l|2003l ) speculated that the 
companion may be evolved and earlier than G5, based on 
aperiodic variability at very low frequencies, similarly to Z 
sources. 

Simultaneous X-ray and optical observations were per- 
formed in 1999 August by the Rossi X-ray Timing Explorer 
(RXTE) and the Sou th African Astron omical Observatory 
(SAAO), respectively ijKong et al.ll2006l ). The orbital period 
was confirmed in the optical observations, but no corre- 
sponding X-ray modulation was found. Neither was there 
any significant X-ray /optical correlation in the light curves. 
Several different two-component spectral models consisting 
of a blackbody and a Comptonized component were suc- 
cessfully fitted to the data. It had already been established 
previously that the spectrum can not be fitted well with a 
single-c omponent b lackbody, power law or bremsstrahlung 
mo del (fSchul3ll999l). 

iLevine et al.l l|2006l ) reported a strong increase in the 
orbital modulation of the X-ray flux in the 2-12 keV band 
since 2005 January, based on over ten years of RXTE All- 
Sky Monitor (ASM) observations. The peak-to-peak ampli- 
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tude of the modulation as a fraction of the average intensity 
was at most 6 per cent from the beginning of the RXTE 
mission in 1996 until 2005 January 19, but was about 18 
per cent thereafter to 2005 May 25 and in four intervals be- 
tween the latter date and 2006 June 9. The modulation was 
more or less independent of energy. 

The distance of GX 9+9 is not well established, but 
it is generally considered to be a Galactic bulge object. 
An estimate of 5 kpc i s almost universally used (from 
Chris tian fc Swankl 119971 . based on visual magnitude, X- 
ray flux during 'burst' episodes and an assu med similarity 
with X1735-44), but occasionally so is 7 kpc ijSchaeferll 19901 : 
IVilhu et alJl2007h. 

Recentlv Hvilhu et all l|2007h studied INTErnational 
Gamma-Ray Astrophysics Laboratory (INTEGRAL) and 
RXTE spectra of GX 9+9 from 2003 and 2004 in the frame- 
work of the spreading layer theory, a model for the luminous 
zone formed on a neutron star surface as a continuation of 
the traditional geometrically thin accretion disc. The spectra 
were fitted with a model consisting of two blackbody compo- 
nents, one of which is multicolour and originates in the inner 
accretion disc, while the other is weakly Comptonized and 
comes from the spreading layer. Estimates of the angular 
extent of the spreading layer on the stellar surface seemed 
to agree with the theory, but there was unexpectedly an ap- 
parent increase of the spreading layer temperature at low 
luminosities. 

This paper is a continuation of the work of IVilhu et al.l 

l|2007h . using a much larger data sample. Section 2 discusses 
the spreading layer theory, Section 3 presents the observa- 
tions used in this work and the data reduction methods, Sec- 
tion 4 introduces the spectral model, while Section 5 includes 
the fitting results, Section 6 some discussion, and Section 7 
the final conclusions. 



2 SPREADING LAYERS 

In non-pulsating neutron star (NS) LMXBs, the magnetic 
field of the neutron star is weak enough 10 s G) not to 
affect the accretion flow, and so the accretion disc can extend 
to the neutron star surface. If the star rotates slower than the 
Keplerian velocity of the infalling matter, the difference in 
kinetic energy must be radiated in a boundary layer between 
the accretion disc and the stellar surface. This energy is of 
the same order as that radiated in the accretion disc, but 
comes from a smaller area closer to the star, which should 
be hotter than the disc and thus produce harder radiation. 

The original acc retion disc model of 
Shak ura fc Sunvaevl dl9T3 | ) led to a bou n dary layer 
(BL) model (see e.g. IPopham fc Sunvae 3 (j200lh and 
references therein) that was thin both radially and lat- 
itudinally, in which the turbulent friction between the 
differentially rotating layers of gas was responsible for 
slowing it down to th e stella r rotational velocity. However, 
llnogamov fc Sunvaevl jl99 9) introduced a new approach 
where the principal friction mechanism is the turbulent 
viscosity between the infalling matter spreading on the 
stellar surface and the slowly moving dense matter beneath 
it. After an intermediate bottleneck zone, the latitudinal 
velocity ve exceeds the radial velocity v r , and the plasma 
forms a distinct spreading layer (SL), dynamically sepa- 




Figure 1. The spreading layer geometry. Matter flows from the 
accretion disc onto the neutron star surface and forms the SL, 
defined as the faster-rotating, radiating part of the spreading flow. 
is the boundary angle of this zone, outside which the surface is 
relatively cool and dark, vg is the latitudinal velocity and v v the 
azimuthal (rotational) velocity. 



rate from the accretion disc. Fig. [T] illustrates the basic 
geometry. The width of the spreading layer, defined as the 
area where the matter rotates faster than the underlying 
surface, increases along with the accretion rate M, and 
consequently the luminosity L. The best-fitting power-law 
dependence to these in the numerically integrated model 
had an index of ~ 0.8. At a certain rate M po i c almost the 
entire surface is covered by the SL. The radial thickness 
and its profile also depend on M; at their maximum in the 
disc plane they are in the range of 0.2-2 km. Unlike the 
traditional boundary layer model, where the luminosity 
has a maximum at the equator, the spreading layer model 
turns out to have luminosity maxima near the outer edges 
of the radiating zone, while at the equator there is a local 
minimum. Radiation from below the spreading flow is 
reprocessed into the blackbody-like SL spectrum in the 
optically thick, higher M case and Comptonized at very 
sma ll M and t. 

[Suk imanov fc Poutanenl l|2006h calculated the spectra 
produced by a spreading layer at different luminosities, and 
seen from different inclination angles; the dependence on 
both of these was slight. The SL spectrum was found to 
be close to a diluted blackbody. Both the effective tem- 
perature Ted and the hardness factor / c , with which it is 
multiplied to get the observed colour temperature T c of the 
black body, increased slightly with SL luminosity. The the- 
ory of llnogamov fc Sunvaevt l|l999l ) was developed by taking 
into account chemical compositions other than pure hydro- 
gen and a general relativity correction to the surface gravity. 
In the original theory, energy release was assumed to happen 
mainly in a thin sublayer at the bottom of the SL. Here it 
was shown that the spectra depend very little on the verti- 
cal structure of the SL when the surface density is high (i.e. 
optically thick case). 

A few papers have compared the theory with obser- 
vations, as well as with the predictions of the traditional 
boundary layer theory. Evidence supporting the SL scenario 
has been found, but the superiority of either model is yet to 
be e stablished conclusiv ely. 

IChurch et aD l|2002l ) compared the SL and BL theories 
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with the results of a survey of LMXBs using the Advanced 
Satellite for Cosmology and Astrophysics (ASCA), which 
were originally published by IChurch fc Balucinska-Churchl 
(2001). They found good agreement with the former at low 
luminosities, but at higher luminosities the blackbody emis- 
sion exceeded predictions by a factor of two to four, suggest- 
ing that radial flow d omina ted the e mitting area. 

Gilfanov et al.l l|2003l ) and iRevnivtsev fc Gilfanovl 



(2006) used Fourier frequency resolved spectroscopy to 
study the short term variability of luminous LMXBs. 
They found that the shapes of the more variable (on 
^second-millisecond time scales) hard component were 
blackbody- like and luminosity-independent, with colour 
temperatures of 2.4 ± 0.1 keV. This was co n sistent 
with the SL model. ISuk imanov fc Poutanenl (2006) 
also compared their theoret ical SL spectra to those of 
IRevnivtsev fc Gilfanovl |2006), finding in five cases out of 
six a soft excess below 10 keV. This was interpreted as 
possibly coming from a classical boundary layer, existing 
between the accretion disc and the spreading layer. 



3 OBSERVATIONS AND DATA REDUCTION 

The basic data set for this work was formed by all avail- 
able pointed observa tions of GX 9+9 by the INTE GRAL 
ijWinkler et al.ll2003l ) and RXTE ijBradt et al.lll993h satel- 
lites between 2002 and 2007. 

For INTEGRAL, we used d ata from both Joi nt Eu- 
ropean X-ray Monitors (JEM-X, iLund et al.l 120031 : JEM- 
X2 before 2004 March, JEM-X1 afterwards) and the IN- 
TEG RAL Soft Gamma-Ray Imager (ISGRI, iLebrun et ail 
2003) - the low energy dete ctor of the Imager on Board the 
INTEGRAL Satellite fIBIS lUbertini et al.l |2003l). All data 
were reduced using the Off-line Science Analysis (OS A) ver- 
sion 7 software. To create average spectra, all pointings with 
more than a few hundred seconds of good JEM-X data and 
with offsets below 5° were used. For the individual spectra, 
to ensure sufficient S/N, we limited the offsets to 4°. 

The RXTE instruments utilized were Proportional 
Counter Unit s (PCUs) and 2 of the Proportional Counter 
Array fPCA. Ijahoda et al.lfl996T ). as the rest are off most 
of the time, using all three Xenon layers, and mostly clus- 
ter A, but for a few observations when it was off cluster 
B of the High Energy X-ray Timing Experiment (HEXTE, 
iRothschild et alj|l998T l. PC A spectra were extracted from 
the Standard-2 mode data, using 16 s binning, and HEXTE 
spectra from the Archive mode data, using 32 s binning to 
match the dwell times. 

Observations over several orbital periods were divided 
into continuous single-orbit viewings of the source. The time 
resolution is thus about 3200 seconds. 

The Good Time Interval (GTI) criteria were the follow- 
ing: PCUs and 2 must be on, the offset from target less 
than 0.02°, the elevation above 10° to eliminate earth oc- 
cupations and atmospheric influence, the parameter value 
ELECTRON2 (a measure of background electron contami- 
nation for PCU 2, usually identical for the other PCUs) be- 
low 0.1 and the time since South Atlantic Anomaly passage 
more than 10 minutes (or negative). The GTI extensions in 
the data files themselves were also applied, except for one 
observation where this resulted in no data at all. 



The PCA deadtime correction factors were calculated 
from the corresponding Standard-1 data. The epoch 5 bright 
background model was used to create the background data 
files, from which the PCA background spectra were ex- 
tracted. The two HEXTE background fields produced by 
the beamswitching were reduced together. 

The final set consisted of 42 JEM-X1+ISGRI, 62 JEM- 
X2+ISGRI and 92 PCA+HEXTE spectra (196 in total), 
spanning 2002 May 1 to 2007 July 4. 



4 THE MODEL 

The spectral modelling was do ne with Xsp e c v. 12 . 4.0, using 
a similar model as the one in IVilhu et al l l|2007l ). It essen- 
tially consists of two modified blackbody components, one of 
which represents the accretion disc and the other the spread- 
ing layer, plus interstellar absorption. The composite model 
is defined as const * wabs(diskbb + compbb). 

LMXB spectral models such as this, based on 
diskbb + compbb, have been dubbed eastern, in contrast to 
the western models based on a basic blackbody from the 
boundary layer and a Co mptonized continu um from a disc 
corona, e.g. bb + comptt |White et al.l ll986 r ): for a compar- 
ison see IPaizis et al.l l|2005l h 

Compbb is a Co mptonized blackbody model after 
iNishimura et al.l ljl986T l. This model is used for the spread- 
ing layer. There are four model parameters: blackbody tem- 
perature kT, electron temperature of the Comptonizing hot 
plasma kT c , optical depth of the plasma r, and the normal- 
ization Nsl = (R-km/ Dio) 2 , where /Zkm is the source radius 
in km (if the source is spherical, which is not the case here; 
see Section 16.211 and Dio is the distance to the source in 
units of 10 kpc. The Comptonization part is accurate for 
photon energies up to the (non-relativistic) electron tem- 
perature, i.e. E <kT c < m c c 2 , and particularly important 
when the optical depth is small (r<3). However, here we 
set kT c to equal the blackbody temperature kT, as we have 
found the Comptonization hard to constrain and weak even 
at best, and assum e it to take place on or near the SL surface. 
IVilhu et all l|2007r i found their best-fitting optical thickness 
to be 0.53, using kT c = 10 keV. 

Diskbb is a multi-temperat ure blackbody model f or ac- 
cretion discs, discussed in e.g. iMitsuda et al.l ll 19841 ) and 
iMakishima et al.l l|l986l ). The model parameters cirG -tin. the 
temperature at the inner disc radius, and the normalization 
-/VDisc = (-Rin/Dio) 2 cos i, where R ln is an 'apparent' inner 
disc radius in km, Dio is the distance as in compbb, and % is 
the inclination of the disc normal from the line of sight. 

The true inner radius of the disc R, which we also 
assume to be the neutron star and spreading layer ra- 
dius (ignoring the thickness of the layer), relates to R\ n as 
R = Rin^n 2 , where £ = (3/7) 1/2 (6/7) 3 w 0.41 is a fac- 
tor that arises from a relativistic boundary condition, which 
causes the disc te mperature to reach its maximum outside 
the inner radius l|Kubota et al.|[l998h . and k is a spectral 
hardness factor, which can be assumed to stay at a bout 1.7 
ijMakishima et al.ll200ol : IShimura fc Takaharalll995l ). There- 
fore setting Rm to 10 km roughly corresponds to an actual 
neutron star radius of 12 km when taking this factor into 
account. 

Wabs is an interstellar photoelectric absorption model 
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using Wisconsi n cross sections llMorrison fc McCammonl 
Il983h and the lAnders fc Ebiharal [jl982f ) relative element 
abundances. The neutral Hydrogen column density nn was 
froze n to 0.196 ■ 10 22 cm~ 2 , calculated from the LAB sur- 
vey ijKalberla et al.|[2005h weighted average of seven points 
within 1° of GX 9+9, by the HEASARC n H web tool. 

The model was also multiplied by a constant, 
energy-independent factor (const) for instrument cross- 
calibrations. It was frozen to the deadtime correction fac- 
tor for the PCA spectra (1.002-1.008), or to unity for the 
JEM-X spectra, and left free for the other instruments. 



Table 1. Best-fit parameters of the average spectra. T; n and 
kT are the inner accretion disc and spreading layer tempera- 
tures, respectively, while A§l is the spreading layer normaliza- 
tion. The accretion disc normalization A^isc = 36l 2 was tied 
to the PCA parameter for all instruments. The Comptonizing 
electron temperature kT c was tied to the SL blackbody temper- 
ature. PCA and JEM-X spectra were fitted together with corre- 
sponding HEXTE and ISGRI spectra, respectively. The bolomet- 
ric model (without absorption) flux was 9.7- 10 -9 erg cm" 2 s _1 , 
2.3- 10 — 9 erg cm -2 s _1 of which was from the compbb component. 
x 2 /dof = 256.44/277 ss 0.93. 
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The energy bands taken into account were 3-27 keV for 
PCA, 4-18 keV for JEM-X1 and JEM-X2, 10-27 keV for 
HEXTE and 12-42 keV for ISGRI; the upper limits were set 
where the background level typically exceeded the signal for 
each instrument. 3.0 per cent systematic errors were used 
for the individual JEM-X+ISGRI spectra, and 2.0 per cent 
errors for the individual PCA+HEXTE spectra. The fitting 
and the la confidence interval calculations were done by the 
Levenberg-Marquardt method. 

The averaged spectra from all the instruments, fitted 
simultaneously, are shown in Fig. [2] In the fitting, the Nuisc 
parameters for the HEXTE, ISGRI and JEM-X datasets 
were tied to the PCA Afoisc parameter. This is because we 
assume the accretion disc to constantly reach the neutron 
star surface, feeding the optically thick spreading layer and 
enabling the blackbody-like hard spectral component to ex- 
ist. This leads to A/bisc being constant, as defined in the 
previous section. The PCA data were used to establish the 
value, as they had the tightest statistical error limits. 

All other parameters for the HEXTE and ISGRI 
datasets were tied to their equivalents in the PCA and JEM- 
X datasets respectively. 

As the JEM-X1+ISGRI and JEM-X2+ISGRI obser- 
vations are from different times than each other and the 
PCA+HEXTE observations, the time-variable parameters 
Ti n , kT (= kT c ), T, and Asl for these spectral groups were 
left independent of each other. 2.5 per cent systematic er- 
rors were used for all the average spectra. The best-fitting 
parameter values are given in Tabled] As can be seen here, 
the JEM-X+ISGRI spectra favoured rather different values 
than the PCA+HEXTE spectra. The unfolded model fitted 
to the average spectra can be seen in Fig. [3] 

The individual spectra were then fitted with the same 
model as for the averaged spectra, but in this case the Nrnac 
parameter was frozen to the best-fitting value obtained from 
the averaged spectra. Taking into account the confidence 
limits (36^2)! the fits were thus done with A/bisc = 34, 36 
and 37. The final confidence limits for the other parameters 
were obtained by summing in quadrature their la errors in 
the Abisc = 36 fits with the differences in the best-fitting 
values between the fits using A?Disc = 36 and those using 
Abisc — 34 and Abisc = 37. 

The Comptonization optical thickness was hard to 
constrain for the individual spectra, so as an initial hy- 
pothesis it was frozen to zero, as suggested by the aver- 
age PCA+HEXTE spectra (which were more constraining 
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Figure 2. Average spectra for the instruments, folded 
const * wabs(diskbb + compbb) model, and data to model ratios. 
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Figure 3. Unfolded model corresponding to Fig. [2] The lower- 
energy family of component curves represents the accretion disc 
component (diskbb) and higher-energy family the spreading layer 
component (compbb). 
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Figure 4. (a) A colour-colour diagram of GX 9+9. The grey area is adapted from lGladstone et al . (2007) and represents all the atoll 
source data; (b) the boxed area of (a); (c) A corresponding colour-luminosity diagram, assuming a distance of 10 kpc and a mass of 1.4 
Mq; (d) the boxed area of (c). 



than the JEM-X+ISGRI spectra). This effectively rendered 
compbb into a basic area-normalized blackbody (bbodyrad). 
Good fits were obtained for nearly all the individual spectra. 
For those with x 2 I dof > 1.2, 19 out of 195, another fit was 
made with r as a free parameter, to see if it would increase 
and if the fit would improve. This was not the case for any 
of the 8 PCA+HEXTE spectra. 5 out of 11 JEM-X+ISGRI 
spectra gained an increase in r, but even at best x 2 1 dof 
only went from 139.17/114 to 132.97/113. It was concluded 
that partial Comptonization in the SL was not significantly 
detected. 

The results from four occasions of partially overlapping 
INTEGRAL and RXTE observations from 2003 October (11 
spectra) were seen to be consistent within error. 

The means for the free parameters in the individual fits 
were T in = 1.787 ± 0.005 keV, kT = 2.50 ± 0.02 keV and 
N s -l = 5.5 ± 0.2; mean x 2 /dof = 0.957 ± 0.014. Sample 
tables of the best-fitting parameters (full tables in the on- 
line version of this paper) can be found in Appendix [B] 



5.2 Fluxes, colours and luminosities 

Response-independent, unabsorbed fluxes were integrated 
over the model spectra over four energy bands: band A = 



3-4 keV, band B = 4-6.4 keV, band C = 6.4-9.7 keV and 
band D = 9.7-16 keV. Hardness ratios were calculated from 
the fluxes, defining th e soft colour as Fb/Ja a nd the hard 
colour as Fp / Fg , as in Done & Gier ihski (2003), a nd s ubse- 
quentlv lGiadstone et al.l ||2007m and lVilhu et ail l|2007l ). The 
colour-colour diagram of Figs. Ufa) and[4jb) shows GX 9+9 
consistently occupying the banana state, as expected. The 
typical correlation seen in GX 9+9 is more evident in the 
PCA+HEXTE results, which had considerably less statisti- 
cal variance. The earliest PCA+HEXTE observations from 
2002 May 1, May 2, June 6 and June 11 form their own 
somewhat distinct track (the purple open squares). 

Bolometric fluxes and compbb (SL) component fluxes 
were calculated, the latter by setting the diskbb normal- 
ization iVbisc to zero. The diskbb fluxes were derived by 
subtracting the compbb fluxes from the corresponding total 
fluxes. 

To calculate bolometric luminosities from the fluxes, a 
value has to be assumed for the distance. Let us consider the 
parameter AT Disc = (7? in /Di ) 2 cos i = 36±2- Using R m = 10 
km, as discussed previously, and the most generally used 
value for the distance, 5 kpc, would result in an inclination 
i of 84.8°. This is unrealistic, as even though t he orbital X- 
ray m odulation has increased since early 2 005 llLevine et a~l 
2006), GX 9+9 is not an eclipsing binary. IPaczvnskil l|l97ll ) 
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relates the radius of the secondary Roche lobe r\ and the 
orbital separation A to the masses of the secondary and the 
neutron star, M\ and M2: 

% = 0.46224 V3 
A \M 1 +M 2 J 

The lack of eclipses constrains this to be < l/(tan i), so 
for a 1.4 Mq neutron star and a 0.2 -0.45 secon dary 
the inclination upper limit is 74°-77° . ISchaefd l|l990h fur- 
ther reasons that the lack of significant dips in the light 
curve implies an inclination less than ~ 70°. On the other 
hand, distances larger than 11 kpc are increasingly unlikely, 
as the binary would be further and further away from the 
Galactic bulge, and its luminosities increasingly above the 
Eddington limit L^dd- Therefore for the rest of this paper, 
we assume a distance of 10 kpc. While being twice the usu- 
ally quoted value, it has the merit of leaving the inclination 
just below 70° without leading to supercritical luminosities, 
and is quite reasonable for a Galactic bulge object. Fig. [5] 
shows the parameter dependences for the best-fitting value 
and confidence limits, as well as the possible distances along 
the line of sight in Galactic context. This new distance es- 
timate does put GX 9+9 in a less populated area in the 
colour-luminosity diagram of Figs. Sic) and[4jd), in the ap- 
proximate range of 0.5-0.9 I/Edd- Fig. [6] shows the relative 
luminosities of the accretion disc and spreading layer com- 
ponents. There is no clear correlation between these, ex- 
cept for the earliest PCA+HEXTE observations from 2002 
May/June, where the accretion disc luminosity was rather 
constant at a little over 0.4 I/Edd- The correlation is smeared 
out by variations in the total luminosity, in a smaller part 
due to orbital modulation. 



50 













* IW '** , *"*** JH:w . R = 16 


km 




***** 










R = 12 km 7 ^ 
























R = 8 km **J;X 






f (a) 


•A '-. 
\\ ■. 

'". 

■■v 


\ \ 
\\ \ 

'■A *'• 

\\ \ 
, '. \ , 





4 6 8 10 12 

Distance [kpc] 




Figure 5. (a) Inclination versus distance for iVbisc = 
(i?i n /Dio) 2 cos i = 36^2 a * three different values of R fa 1.19Ri n 
(curves). The horizontal lines are the maximum inclinations in 
the case of a 0.2 or 0.45 Mq secondary, 77° and 74°, while the 
vertical dashed line denotes the distance of ~ 11 kpc, where the 
high end of the luminosity distribution calculated from the model 
fluxes surpasses the Eddington limit. The region of acceptability 
is shaded; (b) Suggested locations of GX 9+9 along the line of 
sight (white dots) and a cross section of the Galactic environ- 
ment (after Buser 2000). GX 9+9 most likely belongs to either 
the bulge or the thick disc population. 



6 DISCUSSION 



Of all the available spreading layer variables, luminosity 
should be the one most closely correlated with mass ac- 
cretion rate, the driving force behind observed changes not 
due to orbital modulation. Therefore the other spreading 
layer parameters are presented in relation to it. As can be 
seen from Fig. kT and Asl vary considerably within a 
few hours or days, but are close to constant on longer time 
scales, with perhaps a slight indication of a rising trend in 
the former and a decreasing trend in the latter. 



6.1 Spreading layer temperature 

In the models of ISuleimanov fc Poutanenl |2006l ). the effec- 
tive SL temperature T e s decreases slightly with decreasing 
SL luminosity. The observed colour temperature (kT = T c = 
/cTeff) res ults in Fig. ETa) sh ow an opposite trend, as they 
also did in IVilhu et alJ|2007l ). This may be due to an actual 
increase of the effective temperature from some low accre- 
tion rate factor that has not been properly considered in the 
theory, or the luminosity dependence of the hardness factor 
fc = T c /T c g. In the latter case, / c should have a similar 
or somewhat stronge r low-luminos i ty inc rease. The analyti- 
cal formula for / c of Pavlov et al.l ljl99ll ) that was used by 
ISuleimanov fc Poutanenl l|2006t ) instead gives a strong de- 
crease at luminosities below 0.1 I/Edd: 



X 

x .Prk^n XX 




JEM-X1 +ISGRI + 
JEM-X2+ISGRI x 
PCA+HEXTE □ 
PCA + HEXTE, □ 
2002 Moy/June 



X 



□ _ 



0.05 0.10 0.15 0.20 0.25 0.30 0.35 
Spreading layer luminosity [ L SL / L Edd ] 

Figure 6. Luminosities of the two model components. 



3 + 5X\- 3 ft\ 
L 2 o (i) 



where X is the hydrogen mass fraction and the luminos- 
ity L is expressed in Eddington units. This formula has 
been successful in describing ne ar-Eddington lumin osity 
(> 0.9 I/Edd) X-ray burst spectra ijPavlov et al.lll99"lh . but 
may be incorrect in the low-luminosity domain where it has 
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Figure 7. Best-fit spreading layer parameter values vs. time 
(MJD). (a) kT; (b) N SIj . 



not been extensively tested. As can be seen in Fig. E^b), 
where the effective temperature as given by Equation is 
shown, this dependency for the hardness factor only serves 
to highlight the discrepancy between the results and the 
theory, leading to a clear supercritical temperature increase 
towards the lower luminosities. The critical Eddington effec- 
tive temperature TEdd is determined by the balance between 
the surface gravity and the radiative acceleration: 



GM 



R?y/1-Rs/R 



GMc 



asBCTe-ffV 1 - 2GM/c 2 R 



where a c = 0.02(1 + X) m 2 kg -1 is the electron scattering 
opacity, X is the hydroge n mass fraction and o~sb is th e 
Stefan-Boltzmann constant ISuleimanov fc Poutaner] j2006). 
For solar composition plasma (X = 0.7), M = 1.4 Mq and 
R = 12000 m, T Edd ~ 1.92 keV. 

6.2 Spreading layer boundary angle 

As mentioned in Section |4l the spreading layer normaliza- 
tion A^sl is proportional to i?^ m and to the inverse of D\ . 
As the spreading layer is not a spherical source, -Rkm / R. 
To estimate the extent of the SL area to even some degree 



Figure 8. (a) Observed spreading layer colour temperature vs. 
luminosity; (b) effective spreading layer temperature as calculated 
from (a) by dividing with the hardness factor f c (Equation[TJ. The 
RXTE /PCA+HEXTE error bars are mostly within the symbols. 
Also shown are best-fitting power laws (x 2 /dof ~ 2.7) and the 
Eddington temperature for spherical emission. 



of accuracy, we need to relate these with a corrective term; 
this should, in effect, be the ratio between the projection on 
a plane perpendicular to the line of sight of the visible, lu- 
minous SL area Asl and ttR 2 , the corresponding projection 
of the whole neutron star area. 

Here we ignore relativistic light bending, which does 
actually cause the far hemisphere of a neutron star to be 
visible up to an angle of ~20-40° beyond the classical hori- 
zon. Even so, a f ew further s i mplify ing assumptions have to 
be made. As in IVilhu et al.l l|2007t ). the spreading layer is 
approximated by a spherical zone extending from the plane 
of the accretion disc to a certain boundary angle 8; the ac- 
cretion disc is opaque and large enough compared to the 
neutron star to hide the other side of the SL at all inclina- 
tions i < 90° . This approximation probably tends to under- 
estimate the upper boundary angle, as according to theory 
most of the radiation should come from a narrower belt that 
doesn't reach the equator. But as we can only presume to 
see the projected area and not the shape of it, or its position 
on the stellar surface, having two different boundary angles 
as free parameters might give arbitrary results. The mathe- 
matics of the approximation are presented in Appendix [Al 

The results agree well with those of a simple Monte 
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Figure 9. Approximate spreading layer boundary angle vs. lu- 
minosity and best-fitting power law (x 2 /dof 1.5). 

Carlo simulation. For small values of (Nsl, 0), both the ap- 
proximation and the simulation depend only weakly on the 
inclination at i > 60°. The approximation becomes less ac- 
curate the more of the SL on the other hemisphere is visible. 

The resulting boundary angle values are depicted in Fig. 
[9] as a function of spreading layer luminosity. The mean of 9 
was (3.67 ± 0.14)°. The best-fitting power law dependency 
is close to line ar, as opposed to the theor etically predicted 
index of ~ 0.8 lilnogamov fc Sunvae vll99gh . Without proper 
relativistic corrections, the implications on theory remain 
unclear. 



7 CONCLUSIONS 

We have achieved formally successful fits to a large number 
of GX 9+9 spectra with a diskbb + compbb model, inter- 
preting the compbb part as spreading layer emission. The 
results, however, were somewhat unexpected. 

Assuming the model of an accretion disc constantly 
reaching the neutron star surface is valid, GX 9+9 may be 
twice as distant as previously thought, of the order of 10 kpc. 
The distance is suggested by the observed normalization of 
the soft spectral component identified with the accretion 
disc (36! 2 ), together with the upper limit imposed on the 
system inclination (~ 70°) by the lack of eclipses. This dis- 
tance corresponds to a luminosity range of ~ 0.5-0.91/Edd- 
Lack of information on the actual neutron star radius denies 
us accurate boundaries for the distance. 

There was an apparent increase in the temperature of 
the spreading layer component at low spreading layer lu- 
minosities, either due to an effective temperature increase 
from some low accretion rate factor not considered in the 
theory, or incorrect theoretical low-luminosity values for the 
hardness factor f c . 

The spreading layer blackbody parameters varied con- 
siderably on time scales of a few hours or days, but were close 
to constant in the long term, the means being kT w 2.5 keV 
and A r sL ~ 5.5. The Comptonization parameters were found 
to be unnecessary, i.e. no hardening to distort the blackbody 
shape was significantly detected. 

Apart from the low-luminosity temperature rise, the 



results were consistent with spreading layer theory. The 
relatively stable SL colour temperature of ~ 2.5 keV 
is also compa t ible w ith t he 2.4 ±0.1 keV obse r ved i n 
iGilfanov etall (|2003h and iRevnivtsev fc Gilfanovl |2006). 
The aforementioned short-term studies might have missed 
a low-luminosity temperature rise, if it did occur in the 
observed sources. Fourier frequency resolved spectroscopy 
should also be applied to GX 9+9 at various SL luminosities 
to further compare the results of the different approaches. 

There are also some differences in the results of this and 
earlier obser vational studie s on th e subject. Whereas the ob- 
servations in IChurch et all (2002) had good agreement with 
the SL theory at low luminosities and trouble at higher ones, 
in this work the opposite seems to be the case. The observed 
excess below 10 keV reported by ISule" imanov fc Pout anen 
( 2006) was also not detected in our study. 
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Adding the area of the lower ('equatorial') half-ellipse and 
subtracting the area of the upper one, we get an estimate 
of the shaded area in Fig. lAlT b) that should be viable at 
reasonably high inclinations and low spreading layer angles: 



2 



n — 29 sin i 
2^ 



„ / 1 _ . . . . R sin 9 sin i 
2 —R sm 9 sin i 
\2 



tan(# sin i) 



nR . ttR R sin 9 sin i 

H cos i t- r cos 9 cos i 

2 2 tan(0 sin i) 

Now A^sl is solved to single decimal accuracy and compared 
to observations to find 9: 



R /#sinj sin 9 sin i cos i 
D\ Q V 7r 7r tan((9 sin i) 2 

sin 9 sin i 



2 tan(# sin i) 



-cos 9 cos i 



R=12 (km), Dw = 1.0 (= 10 kpc) and i = 69.0 (°) were 
used for the calculations (Fig. [9}. 



APPENDIX A: SL BOUNDARY ANGLE 
APPROXIMATION 

The extreme case closer to the likely inclination of the sys- 
tem (60-70°; see Section T5.2P is the case of i ~ 90°. Viewed 
nearly edge-on, the accretion disc allows us to see only a 
semi-circular area of the stellar surface. To get Asl from 
this, we have to subtract the polar segment, i.e. the area of 
the sector whose central angle is n — 29go, minus the triangle 
within the spreading layer: 

Asl = — ( — — — tt-R -2--R sm 6»9o cos 9O ) 

One way to approximate the projected area at intermedi- 
ate inclinations close to 90° is to start with the equation 
above and add a few corrective terms. As depicted in Fig. 
lAlT a). 'vertical' lines and angles are multiplied by the sine 
of the inclination, while 'horizontal' lines are multiplied by 
the cosine. This doesn't apply to the semimajor axis of the 
half-ellipse that forms the upper boundary of the SL projec- 
tion; it has to be derived from the projected angle and the 
minor cathetus: 

R sin 9 sin i Rsin9sini 

tan(# sm i) = <4> x = 7-7 r- 

x tan(t/ sm 1) 
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Figure Al. The geometry of the projected visible spreading layer area approximation, (a) cross-section side view; (b) observer view. 

APPENDIX B: BEST-FITTING PARAMETERS This paper has been typeset from a TpX/ MpX file prepared 

by the author. 
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Table Bl. Observation ID:s, dates, Modified Julian Dates, best-fitting parameter values, bolometric model fluxes (in 10 9 erg cm 2 
s _1 ) and fit statistics of the RXTE /PCA+HEXTE spectra. 
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oz -0.04 


— 0.5 


9.33 


52 


.54/55 


70022- 


-02- 


-04- 


-01 


01/13/03 


52652 


.55 


1 


7R+0.03 
'"-0.02 


2 


qc+0.03 
••-' o -0.02 


6.5+ n '5 

—0.7 


9.60 


39 


.69/55 


70022- 


-02- 


-04- 


-00 


01/13/03 


52652 


93 


1 


7C-+0.03 
•'°-0.02 


2 


41+O.O2 
-0.02 


— 0.6 


10.17 


43 


.46/55 


70022- 


-02- 


-04- 


00 


01/13/03 


52653 


.00 


1 


77+0.03 
' ' -0.02 


2 


41+O.O3 
•^ 1 -0.02 


—0.6 


10.11 


42 


.61/55 


70022- 


-02- 


-04- 


■00 


01/13/03 


52653 


.06 


1 


7C-+0.03 
' °-0.02 


2 


40+O.O2 
'^"-0.02 


—0.6 


10.43 


36 


.93/55 


70022- 


02- 


-04- 


-00 


01/13/03 


52653 


.12 


1 


70+O.O3 
'•-'-0.02 


2 


49+0.02 

•^ z -o.oi 


9.6+ n '4 

—0.6 


10.62 


50 


.57/55 


70022- 


-02- 


-04- 


-02 


01/14/03 


52653 


.23 


1 


7K+0.03 
'"-0.02 


2 


qO+0.02 

°°-0.02 


7.8 +0 ,i 

' ^—0.6 


10.16 


39 


.27/55 


70022- 


02- 


■04- 


-02 


01/14/03 


52653 


.29 


1 


7C-+0.03 
,o -0.02 


2 


41+O.O2 
'^ 1 -0.02 


8.3+S1 

—0.6 


10.35 


32 


.18/55 


70022- 


-02- 


■04- 


-02 


01/14/03 


52653 


.35 


1 


7^+0.03 
'"-0.02 


2 


40+O.03 
•^"-0.02 


6.9+°i 

—0.6 


9.95 


54 


.71/55 


70022- 


-02- 


■07- 


-00 


09/16/03 


52898 


.92 


1 


7Q+0.03 
■ ,y -0.02 


2 


qo + 0.03 

■"-0.02 


1 ^—0.7 


10.67 


59 


.59/55 


70022- 


-02- 


■07- 


-01 


09/16/03 


52898 


.97 


1 


7C-+0.03 
,o -0.02 


2 


49+0.02 
'^ z -0.01 


W.2+°i 
—0.6 


11.13 


36 


.44/55 


70022- 


-02- 


-07- 


02 


09/17/03 


52899 


.03 


1 


72+0.03 
' S -0.02 


2 


49+0.02 
-0.02 


"—0.6 


11.08 


47 


.17/55 


ouuzu- 


-U 7 L- 


01- 


-01 


10/09/03 




.0 i 


1 
1 


oo+0. 03 
oz -0.02 


A 


40+O.O5 
' 48 -0.03 


J — 0.6 


10.21 


52 


.36/55 


80020- 


-04- 


-01- 


■00 


10/09/03 


52921 


.66 


1 


09+O.O3 
-0.02 


2 


44+0.05 
'^-0.04 


40+0.5 


10.10 


42 


.59/55 


80020- 


-04- 


-02- 


00 


10/11/03 


52923 


.09 


1 


09+O.O3 
oz -0.02 


2 


46+0-04 

*"-0.03 


44+0.4 
-0.6 


10.32 


45 


.71/55 


80020- 


-04- 


-02- 


■01 


10/11/03 


52923 


.16 


1 


01 +0.03 
8i -0.01 


2 


^9+0.05 
• oz -0.03 


q 9+O.3 
°- z -0.5 


9.72 


64 


.07/55 


80020- 


■04- 


-03- 


-00 


10/13/03 


52925 


.54 


1 


oo+0. 03 
oo -0.01 


2 


crp+0.06 
OD -0.04 


2 q+03 
z ' y -0.5 


10.19 


55 


.49/55 


80020- 


-04- 


-04- 


-00 


10/13/03 


52925 


.91 


1 


oo+0. 03 
oo -0.02 


2 


44+0.03 
'^-0.03 


6 3+ ' 5 

D.O_ 7 


11.29 


48 


.32/55 


70022- 


-02- 


-08- 


-00 


01/16/04 


53020 


.05 


1 


01 +0.03 

-0.02 


2 


4^+0.04 
•^°-0.03 


4 8+ ' 4 


10.23 


58 


.02/55 


70022- 


-02- 


■08- 


00 


01/16/04 


53020 


.11 


1 


oo+0. 03 
8U -0.02 


2 


49+0.03 
'^ z -0.02 


9+O.5 


11.30 


52 


.09/55 


70022- 


-02- 


■08- 


-02 


01/16/04 


53020 


.18 


1 


09+O.O3 
oz -0.02 


2 


47+0.04 
*'-0.03 


4 4+0-4 


10.44 


36 


.49/55 


70022- 


-02- 


■07- 


■04 


01/16/04 


53020 


.90 


1 


7O+0.03 
'"-0.02 


2 


40+O.O2 
'^°-0.02 


U 7+0-5 
L1 - '-0.7 


12.42 


42 


.57/55 


70022- 


-02- 


■07- 


■03 


01/17/04 


53021 


.03 


1 


01 +0.03 
-0.02 


2 


41+O.O3 
•^ 1 -0.02 


6 4 +0 ' 5 

"■ 4 -0.7 


10.74 


36 


.74/55 


70022- 


-02- 


-07- 


-03 


01/17/04 


53021. 


.10 


1 


70+0. 03 
' -0.02 


2 


qo+0.03 
- 3 °-0.02 


74+0.5 
'-0.7 


10.66 


44 


.78/55 


70022- 


-02- 


■07- 


■03 


01/17/04 


53021 


.16 


1 


01 +0.03 
ol -0.02 


2 


44+0.03 
*^-0.02 


5 6+0-4 
°- b -0.6 


10.66 


41 


.42/55 


70022- 


02- 


-07- 


■05 


01/21/04 


53025 


.95 


1 


09+O.O3 
-0.02 


2 


49+0.04 
'^ z -0.03 


c 7+0.5 
°- '-0.7 


10.71 


40 


.98/55 
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Table Bl - continued Observation ID:s, dates, Modified Julian Dates, best-fitting parameter values, bolometric model fluxes (in 10 9 
erg cm~ 2 s _1 ) and fit statistics of the RXTE /PCA+UEXTE spectra. 



Obs. ID 






MM/DD/YY 


MJD 


Tin [keV] 


kT [keV] 




N S -L 




X 2 /dof 


70022- 


02- 


■08- 


-01 


01/22/04 


53026 


02 


1 


on+0.03 
SU -0.01 


2 


46 +0.04 
*°-0.03 


4 


o+°i 

— 0.5 


9.85 


58.17/55 


70022- 


-02- 


■08- 


■01 


01/22/04 


53026. 


.09 


1 


oi +0.03 
•° -0.02 


2 


t-Q+0.05 
OU -0.03 


3 


O+0.4 
-0.6 


9.99 


47.64/55 


80020- 


-04- 


-05- 


-00 


02/18/04 


53053. 


36 


1 


oc+0.03 
'°°-0.02 


2 


n +0.05 

d -0.04 


3 


Q+0.4 

'^—0.6 


10.87 


42.57/55 


80020- 


-04- 


■05- 


■00 


02/18/04 


53053. 


.43 


1 


oi +0.03 
° -0.02 


2. 


49+0.04 
'^-0.03 


4 


4+0-4 

—0.6 


10.05 


53.55/55 


80020- 


-04- 


-06- 


-00 


02/20/04 


53055 


63 


1 


oi +0.03 

'° -0.02 


2 


41 +0.03 
^ -0.02 


6 


O+0.5 

• y -0.7 


10.92 


43.25/55 


70022- 


-02- 


09- 


-03 


04/26/04 


53121. 


.88 


1 


7O+0.03 
■ ' -0.01 


2 


4 7+0.04 
'*'-0.03 


3 


7+0.3 

' —0.5 


9.51 


76.81/55 


70022- 


-02- 


-09- 


■03 


04/26/04 


53121. 


.94 


1 


oi +0.03 
Si -0.02 


2 


A q+0.04 
■*°-0.03 


4 


q+0-4 

' — 0.6 


10.22 


44.76/55 


70022- 


-02- 


■09- 


■04 


04/27/04 


53122. 


.01 


1 


a n+0.03 
SU -0.02 


2 


47+0.02 
-0.02 


8 


9+0.4 

— 0.6 


11.52 


53.88/55 


70022- 


-02- 


09- 


■02 


04/27/04 


53122. 


.82 


1 


on+o.03 

SU -0.01 


2 


-9+0.05 
oz -0.03 


3 


4+0.3 

-0.5 


9.66 


71.75/55 


70022- 


-02- 


09- 


-01 


04/27/04 


53122. 


.93 


1 


fi n+0.03 
SU -0.01 


2 


c-n+o.04 

ou -0.03 


3 


6+0-3 

"—0.5 


9.76 


75.82/55 


70022- 


-02- 


09- 


■02 


04/27/04 


53122. 


99 


1 


oi +0.03 
•° -0.02 


2 


4 7+0.04 
'^'-0.03 


4 


q+0.4 
— 0.6 


10.16 


52.45/55 


70022- 


-02- 


-09- 


■00 


04/29/04 


53124. 


70 


1 


S n+o.03 

SU -0.02 


2 


45+0.04 


4 


5+„^ 
-0.6 


10.00 


43.28/55 


70022- 


-02- 


09- 


-00 


04/29/04 


53124. 


.77 


1 


q-, +0.03 
Si -0.02 


2 


44 +0.04 
* -0.02 


5 


1+0.4 

^-o.e 


10.36 


54.08/55 


70022- 


02- 


-09- 


■00 


04/29/04 


53124. 


83 


1 


7Q+0.03 

•< y -o.oi 


2 


4B+0.04 

n +0.05 
J± -0.03 


3 


7 +0.3 
" ' —0.5 


9.58 


66.55/55 


70022- 


02- 


09- 


00 


04/29/04 


53124. 


.89 


1 


on+0.03 
SU -0.01 


2 


3 


5+0.3 

-0.5 


9.71 


57.94/55 


70022- 


-02- 


09- 


■00 


04/29/04 


53124. 


.95 


1 


on+0.03 
• 8U -0.01 


2 


47+0.04 
'*'-0.03 


4 


i+0.4 

— 0.5 


9.88 


45.10/55 


70022- 


-02- 


-10- 


-00 


06/24/04 


53180. 


69 


1 


QQ + 0.03 

•° -0.02 


2 


cn+o.04 

ou -0.03 


3 


7 +0.3 
' —0.5 

■6+°-5 

— 0.7 


10.09 


43.67/55 


70022- 


02- 


■09- 


-05 


06/24/04 


53180. 


.78 


1 


vq+0.03 
■' S -0.02 


2 


46+0-02 
™-0.02 


9 


11.71 


46.63/55 


80020- 


-04- 


■07- 


■00 


08/26/04 


53243. 


70 


1 


oc+0.03 
•°°-0.02 


2 


ci +0.04 
oi -0.03 


5 


9+0.5 
— 0.7 


11.45 


42.32/55 


80020- 


-04- 


■07- 


-01 


08/26/04 


53243. 


.77 


1 


QQ+0.03 

•°°-0.02 


2 


A c-+0.03 
/io -0.02 


8 


q+0.6 
• — 0.8 


12.10 


37.52/55 


80020- 


-04- 


■07- 


■02 


08/26/04 


53243. 


.84 


1 


QQ + 0.03 

°°-0.02 


2 


A q+0.03 
^°-0.03 


7 


9+O.6 

— 0.8 


11.56 


57.37/55 


80020- 


-04- 


■07- 


-03 


08/26/04 


53243. 


91 


1 


Qc+0.03 
'°°-0.02 


2 


4Q+0.05 
■^°-0.03 


4 


.4+0-5 

— 0.7 


10.95 


41.51/55 


80020- 


-04- 


-08- 


■00 


09/04/04 


53252. 


.69 


1 


oq+0.03 
oo -0.01 


2 


ri+0.06 
u -0.04 


2 


Q+0.3 

— 0.4 


10.03 


107.10/55 


80020- 


-04- 


-08- 


■02 


09/04/04 


53252. 


.75 


1 


Q4+0.03 
' S4 -0.02 


2 


c 4 +0.06 
d ^-0.04 


3 


q+0.3 
u -0.5 


10.29 


71.34/55 


80020- 


■04- 


■08- 


■01 


09/05/04 


53253. 


53 


1 


QQ + 0.03 

■ SJ -0.02 


2 


45+0.03 

-0.02 


6 


6+0.5 
•"-0.7 

0+O.5 
• S -0.7 


11.32 


61.01/55 


80020- 


-04- 


-08- 


-01 


09/05/04 


53253. 


60 


1 


QQ + 0.03 

•° z -0.02 


2 


46+0- 03 
+ D -0.02 


6 


11.34 


32.04/55 


80020- 


-04- 


-09- 


-00 


09/11/04 


53259. 


.98 


1 


QK+0.03 
' Sb -0.02 


2 


47+0.04 

'*'-0.03 
r;c;+0.05 
OO -0.03 


5 


0+0.5 


11.33 


54.09/55 


92415- 


-01- 


01- 


-02 


06/20/06 


53906. 


.35 


1 


77 +0.03 
■' '-0.01 


2 


2 


O+0.3 
• 8 -0.4 


8.93 


75.65/54 


92415- 


-01- 


-01- 


-02 


06/20/06 


53906. 


.41 


1 


7K+0.03 
' D -0.01 


2 


-9+0.04 
oz -0.03 


3 


i+0.3 

• 1 -0.4 


8.78 


53.97/54 


92415- 


-01- 


-02- 


■00 


09/01/06 


53979. 


38 


1 


S9+0. 03 
•° -0.02 


2 


4Q+0.04 
■ 48 -0.02 


5 


9+0.4 
' z -0.6 


10.74 


63.69/55 


92415- 


-01- 


-02- 


■00 


09/01/06 


53979. 


44 


1 


79+0.03 
' 8 -0.01 


2 


/ic+0.03 
^ J -0.02 


5 


O+0.4 
• 8 -0.5 


10.11 


64.51/55 


92415- 


-01- 


02- 


-00 


09/01/06 


53979. 


51 


1 


7Q+0.03 
• ,y -0.02 


2 


n +"."2 
dl -0.02 


9 


4+0-4 
*-0.6 


12.11 


60.49/55 


92415- 


■01- 


02- 


-00 


09/01/06 


53979. 


.57 


1 


QQ + 0.03 

•°°-0.02 


2 


c-n+o- 04 

ou -0.03 


4 


4+0-4 

• 4 -0.5 


10.60 


51.11/55 


92415- 


■01- 


02- 


-00 


09/01/06 


53979. 


.63 


1 


77 +0.03 
■' '-0.01 


2 


44+0.03 
-0.02 


5 


+0.4 


9.66 


56.65/55 


93406- 


09- 


01- 


-00 


07/04/07 


54285. 


.73 


1 


o 7+0. 03 
' 8 '-0.01 


2. 


74+0.07 
'*-0.05 


1 


7+0.2 

'-0.3 


10.66 


89.47/55 



Table B2. Science Window ID:s, dates, Modified Julian Dates, best-fitting parameter values, bolometric model fluxes (in 10 9 erg cm 2 
s _1 ) and fit statistics of the INTEGRAL/ JEM-X2+ISGRI spectra. 



ScW ID MM/DD/YY MJD T in [keV] kT [keV] N SL F Bo i X 2 /dof 



5600760010 


03/31/03 


52729. 


.80 


1 


7CT+0.05 
■'°-0.04 


2 


34+0; 


10 

09 


13 - 6+ 3 


11.72 


126. 


25/114 


5600770010 


03/31/03 


52729. 


.82 


1 


77+0.05 
■' -0.04 


2 


361°; 


17 
14 


8 4+ 3 -* 

° -2.8 


10.47 


109. 


.34/114 


5601000010 


04/01/03 


52730. 


33 


1 


QQ + 0.04 

'° -0.03 


3 


.081°; 


65 
47 


1 0+ 13 


10.82 


101 


.08/114 


5601010010 


04/01/03 


52730. 


35 


1 


09+O.O4 
•° -0.02 


2 


.691° 


18 
16 


3 6+ 14 
— l.i 


10.65 


101 


.08/114 


5601020010 


04/01/03 


52730. 


.38 


1 


01 +0.05 
■ Si -0.05 


2 


.421°; 


14 
12 


11 9+ 4 ' 3 
ll-»_3.4 


12.68 


114. 


.03/114 


5601090010 


04/01/03 


52730. 


.54 


1 


04+0.02 
■°*-0.03 


3 


.78+1 


86 
09 


l+°- 5 


9.11 


124. 


.69/114 


5900290010 


04/08/03 


52737. 


.81 


1 


00+O.O2 
■ Sb -0.02 


3 


.99+1 


28 
25 


2+°- 5 
u ' z -0.2 


9.67 


128. 


63/114 


5900300010 


04/08/03 


52737. 


.84 


1 


oo+0. 03 
°°-0.03 


2 


.78^ 


73 
59 


u -'-0.6 

8 9+ 15 ' 1 


10.12 


108. 


.84/114 


5900310010 


04/08/03 


52737. 


86 


1 


77+O.O6 
■' '-0.10 


2 


00l«; 


43 
23 


9.17 


115. 


37/114 


6100920010 


04/16/03 


52745. 


.26 


1 


7n +0.06 
'' u -0.07 


2 


.20+.0; 


18 
.15 


13 9+ 7 ' 7 


10.13 


124. 


72/114 


6200730010 


04/18/03 


52748. 


.00 


1 


74+O.O6 
•' 4 -0.07 


2 


291° 


33 
24 


6-9±l;l 


9.13 


126. 


12/114 


6200740010 


04/19/03 


52748. 


.02 


1 


77+0.05 
■' '-0.05 


2 


.1310 


14 
12 


H.5l!1 


10.23 


119. 


.04/114 


6200750010 


04/19/03 


52748. 


.04 


1 


69+0-06 
• -0.06 


2 


321"' 


15 
13 


12 -6l^'. 8 


10.31 


120. 


.18/114 
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Table B2 - continued Science Window ID:s, dates, Modified Julian Dates, best-fitting parameter values, bolometric model fluxes (in 
1CT 9 erg cm" 2 s" 1 ) and fit statistics of the INTEGRAL /JEM-X2+ISGRI spectra. 



ScW ID 


MM/DD/YY 


MJD 


T in [keV] 


kT [keV] 


iV S L 


TBol 


X 2 /dof 


6200890010 


04/19/03 


52748. 


.35 


1 


--.q + 0.05 

• '°-0.05 


2 


■08±g;i* 


11.7+5-* 

— 4.2 


9.41 


98.48/114 


6200900010 


04/19/03 


52748 


.37 


1 


7C+0-04 
'°-0.03 


2 




8.3+ 2 ^ 

— 2.3 


9.56 


117.49/114 


6200910010 


04/19/03 


52748. 


.41 


1 


q-i +0.05 
• Oi -0.04 


2 


o 4 +0.24 


6.2 +4 . p 

— 2.9 


10.28 


113.41/114 


6200920010 


04/19/03 


52748. 


43 


1 


70+ 0.06 

■ '^-o.oe 


2. 


Qq + 0.20 

zo -o.i6 


10.3+5-2 

— 4.0 


9.78 


125.70/114 


6300330010 


04/20/03 


52749. 


.83 


1 


7C+0-04 
/O -0.04 


2. 


qo+0.25 
-0.20 


5.6+" 

— 2.6 


9.02 


108.19/114 


6300340010 


04/20/03 


52749 


.85 


1 


70 + 0. 04 
• '^-O^ 


2. 


iq+0.19 
iy -0.16 


— 2.9 


8.65 


115.41/114 


6300350010 


04/20/03 


52749. 


,87 


1 


7^+0.04 
•' D -0.03 


2. 


cc+0.24 
OD -0.19 


3.5+V! 

— 1.4 


9.15 


115.12/114 


6300360010 


04/20/03 


52749 


90 


1 


7Q+0.04 
' ' y -0.03 


2. 


77+0.46 

' '-0.35 


1.6+" 

— 0.9 


8.98 


113.74/114 


6300500010 


04/21/03 


52750. 


.21 


1 


71 +0.04 
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1 
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7.73 
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1 
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04/23/03 
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1 
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9.8+ 12 q 8 
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09 


1 


or+0.04 
' oo -0.04 


2 




5.4+ 3 -t 

— 2.3 


11.07 


95.06/114 


12100320010 


10/11/03 


52923. 


.11 


1 


07+0. 05 
' S '-0.05 


2 




5.5+" 

— 3.3 


11.21 


134.08/114 


12100330010 


10/11/03 


52923. 


.13 


1 


qi +0.03 
yl -0.02 


2 


QC-+0.41 
yo -0.35 


1.3+i"? 


11.41 


106.52/114 


12100410010 


10/11/03 


52923. 


.31 


1 


Qfi+0.03 
° b -0.02 


3 


m +0.36 
u± -0.31 


+ 11 

1*^—0.6 


10.61 


105.37/114 


12200060010 


10/13/03 


52925 


53 


1 


Q 17+0.04 

• 8 '-0.04 


2. 


fi4 +0.55 
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cq+0.27 
do -0.22 
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90 +0.04 
yu -0.03 
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16500640010 
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.60 


1 


97+0.04 
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Table B3. Science Window ID:s, dates, Modified Julian Dates, best-fitting parameter values, bolometric model fluxes (in 10 9 erg cm 2 
s" 1 ) and fit statistics of the INTEGRAL/ JEM-X1+ISGRI spectra. 



ScW ID 


MM/DD/YY 


MJD 


T in [keV] 


kT [keV] 








^Bol 


X 2 /dof 


ZolUUUzUUlU 
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Uy / Uo/ U4 
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oq+0.04 
• oo -0.04 
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iq+0.48 
la -0.43 
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1 
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3 
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1 
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.8 
.6 


1 n 98 
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ZolUU4oUUlU 
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Uy/ U4/U4 


53252 


.82 


1 


oq+0.03 
• 88 -0.02 


2 
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a -0.28 
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1 
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.7 
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lU.yo 
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53252 
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3 
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.70 


1 


79+0.04 
' ,z -0.03 


2 


48+0-13 


5. 
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g-,+0.18 
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5. 




2 
.6 


9.97 
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3 
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1 
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.6 
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.88 


1 
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2. 
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.9 
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1 
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2 
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3 
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